DISTRIBUTION  STATEMENT  A:  Distribution  approved  for  public  release;  distribution  is  unlimited 


Mixing,  Internal  Waves  and  Mesoscale  Dynamics  in  the  East  China  Sea 

Iossif  Lozovatsky  and  Harindra  J.S.  Fernando 

Center  for  Environmental  Fluid  Dynamics,  Department  of  Mechanical  &  Aerospace  Engineering, 
Ira  Fulton  School  of  Engineering  at  Arizona  State  University,  Tempe,  AZ  85287-9809 
*Phone:  480-965-5597,  480-965-2807,  email:  i.lozovatsky@asu.edu  and  J.Fernando@asu.edu 

Grant  Number:  N00014-05-1-0245 


LONG  TERM  GOAL 

The  long-term  goal  of  our  research  program  is  to  better  understand  and  quantify  the  relationships 
between  mesoscale  dynamics,  internal  waves,  turbulence  and  topographic  features  of  shallow,  tidally- 
affected  regions  of  oceans. 

OBJECTIVES 

The  objectives  of  the  project  during  2009  were: 

(i)  to  analyze  the  intermittency  of  turbulence  in  a  tidal  bottom  boundary  layer  (TBBL)  on  a 
shallow  shelf  in  the  western  (Chinese)  sector  of  the  East  China  Sea  (ECS), 

(ii)  to  continue  investigations  on  the  turbulent  kinetic  energy  (TKE)  dissipation  rate  very  near  the 
ocean  floor,  in  particular  the  role  of  small-scale  local  bathymetry  in  small-scale  mixing.  The 
emphasis  was  on  the  area  influenced  by  Changjiang  (Yangtze  River)  Diluted  Water  (CDW).  , 
and 

(iii)  to  study  mixing  efficiency  in  stratified  flows  based  on  data  available  from  atmospheric  and 
oceanic  field  campaigns. 

APPROACH 

Analysis  of  field  data  collected  in  2005  and  2006  in  ECS  during  the  research  cruises  of  the 
Ocean  University  of  China  (OUC)  and  Korea  Ocean  Research  and  Development  Institute  (KORDI). 

WORK  COMPLETED 

Study  of  the  near-bottom  turbulence  in  Chinese  sector  of  ECS:  friction  velocity,  dissipation 
rate,  and  intermittency. 

Study  of  mixing  efficiency  in  stratified  flows. 
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RESULTS 


A.  Turbulence  Intermittency  in  a  Tidal  Flow  near  the  Seafloor  on  a  Shallow 
Shelf  of  the  East  China  Sea 

As  has  already  been  reported  [Lozovatsky  et  al.  2008  a,b],  our  collaborator  and  subcontractor, 
OUC,  conducted  hydrophysical  measurements  in  the  ECS  in  the  close  proximity  (1.2  km)  to  the  shore 
of  the  Jiaozhou  Bay  (A  =  36.04°N,  cp  =  120.32°E).  The  higher-order  (q)  structure  functions  of  vertical 
velocity  fluctuations 

(a  w;)=C,(e,ryM,  (1) 

(transverse  structure  functions  or  TSF)  were  employed  to  study  the  characteristics  of  turbulence 
intermittency  in  a  reversing  tidal  flow  on  a  19-m  deep  shallow  shelf  of  ECS.  Eq.(l)  represents  scaling 
assumption  of  a  refined  similarity  hypothesis  [Kolmogorov  1962;  Obukhov  1962]  that  accounts  for 
fluctuations  of  the  dissipation  rate  s  at  the  scale  r  of  locally  isotropic  turbulence  in  inertial-convective 
and  viscous  dissipative  spectral  subranges.  A  classical  lognormal  single  parameter  (//)  model 

=  4, =  f  “  ~~  3^)  (2) 

for  intermittency  and  a  universal  multifractal  approach  (specifically,  the  two-parameter  (Ci  and  a)  log- 
Levy  model) 

£(?)  =  £>  =f — ^r 
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were  employed  to  analyze  the  TSF  exponent  c(q)  in  tidally-driven  turbulent  boundary-layer  and  the 
parameters  ju,  C\  and  a  were  estimated.  Measurements  from  a  downward-looking  bottom-mounted 
ADV,  positioned  0.45  m  above  the  sea  floor  were  utilized,  and  the  data  spanned  two  semidiurnal  tidal 
cycles..  To  analyze  the  intermittency  in  tidal  boundary  layer  turbulence,  25-hour  records  of  the  ADV 
current  components  were  subdivided  into  44  segments.  Each  segment  contained  2 15  =  32628  individual 
samples  (time  interval  ~  34  min),  which  was  relatively  long  to  provide  sufficient  multiplicative 
averaging  of  higher-order  TSF  calculations  and  to  minimize  errors  to  an  acceptable  level.  The  mean 
TKE  dissipation  rate  snb  at  each  segment  i  was  estimated  from  the  inertial  subrange  of  the  frequency 

counterpart  of  the  Kolmogorov  wave  number  spectra.  Taylor’s  “frozen  turbulence”  hypothesis  was 
employed  to  transform  w(t)  at  each  segment  to  the  spatial  series  w{x).  The  applicability  of  Taylor’s 
hypothesis  was  tested  by  calculating  the  ratio  rmsi(w')/Ul  ;  it  never  exceeded  2.5%  at  segments  close 

to  high  and  low  tides  (minimum  advecting  velocity)  and  mostly  was  below  1%.  Possible  impacts  of 
Taylor  hypothesis  on  TSF  exponents  were  carefully  scrutinized. 
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During  energetic  flooding  tidal  phases  (high  integral  Reynolds  numbers  Ret„i),  the  parameters 
of  intermittency  models  approached  the  mean  values  of  p  ~  0.24,  C\  ~  0.15,  and  a  ~  1.5,  which  are 
accepted  as  the  universal  values  for  fully-developed  turbulence  at  high  Re.  With  the  decrease  of 
advection  velocity,  and  therefore  Re,  p  and  C\  increased  up  to  p~  0.5  -  0.6  and  C\  ~  0.25  -  0.35,  but  a 
decreased  to  about  1.4  (Fig.  1). 
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Fig.  1.  The  parameters  of  multifractal  (Eq.  3,  Cj  and  a)  and  log-normal  (Eq.  2,  p)  intermittency  models 
shown  with  95%  confidence  bounds.  The  turbulent  integral  (Re)  Reynolds  number  is  mainly  in  phase  with  a 
and  out  of  phase  with  Cj  and  p  .  Two  periods  of  flooding  currents  with  the  highest  Re,  recorded  are 

highlighted. 

The  dependence  of  p,  C\,  and  a  on  internal  flow  parameters,  such  as  the  microscale  Reynolds 
number  R^w,  which  is  a  variant  of  Rx  =  rms(u')A/v ,  is  presented  in  Fig.  2.  Despite  high  scatter,  the 

data  show  a  general  increase  of  p  and  C\  with  decreasing  Rxw  (Fig.  2  a,b).  The  best  least-square  power 
trends  are  given  with  95%  confidence  bounds  and  coefficients  of  determination  r  .  The  empirical 
functions  so  obtained  are 

p  =  0.23  +  i?^'8  and  C,  =  0. 13  +  R^  ,  (4) 
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Both  estimates  are  close  to  the  universal  values  of  /}  and  Cl  expected  from  turbulence 
intennittency  models  at  high  Reynolds  numbers.  A  weak  (statistically  insignificant)  dependence  of  a 
on  Rzw  is  shown  in  Fig.  2c.  The  mean  value  a  =1.53  ±0.39,  however,  matches  well  with  the 
“universal”  a  =  1.5  - 1.55  cited  by  Seuront  et  al.  [2005]  for  log-Levy  multifractal  intennittency  model. 
It  should  be  noted  that  the  analysis  here  was  based  on  TSF,  which  could  be  compared  with  the 
laboratory  data  of  Hao  et  al.  [2008],  as  shown  in  Fig.  2  by  stars. 

The  results  explain  the  reported  disparities  between  the  smaller  “universal”  values  of 
intennittency  parameters  //  and  C i  (mostly  measured  in  laboratory  and  atmospheric  high-Reynolds 
number  flows)  and  those  (//  =  0.4  -  0.5)  reported  for  oceanic  stratified  turbulence  in  the  pycnocline, 
which  is  associated  with  relatively  low  local  Reynolds  numbers  Rz,v. 

The  scaling  exponents  c(2)  of  the  second  order  TSF,  relative  to  the  third  order  SF,  was  also 
found  to  be  a  decreasing  function  of  Rzw,  approaching  the  classical  value  of  2/3  at  very  high  R^w.  A 
larger  departure  from  the  universal  turbulent  regime  at  lower  Reynolds  numbers  could  be  attributed  to 
the  higher  anisotropy  and  associated  intennittency  of  under-developed  turbulence. 


B.  The  Dissipation  Rate  near  the  Seafloor  on  a  Tidal  Shelf  in  the  East  China  Sea 

The  profiles  of  the  TKE  dissipation  rate  were  measured  in  the  Changjiang  (Yangtze  River) 
Diluted  Water  of  the  East  China  Sea  (CDW  region:  A  =  30°49’N,  cp  =  122°56’E;  the  mean  depth  Hb  = 
38  m)  by  a  microstructure  profiler  (MSS)  during  two  semidiurnal  cycles.  Basic  characteristics  of  tidal 
flow,  stratification,  and  microstructure  in  the  region  were  analyzed  and  reported  by  Lozovatsky  and 
Fernando[ 2008].  In  2009,  we  continued  investigations  of  CDW  regional  dynamics,  focusing  on  near¬ 
bottom  turbulence.  The  time  variation  of  the  velocity  vector  (its  magnitude  U  and  direction  (f>)  at  Q  = 
0.3  mab  are  shown  in  Fig.  3,  where  Q is  a  distance  in  meters  above  the  bottom. 


Fig.  3.  The  near-bottom  (g=  0.3  mab )  temperature  T„b,  direction  <j>,  and  magnitude  U  of  the  mean  current. 
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The  magnitude  of  the  velocity  vector  was  affected  by  oscillations  having  a  quarter  diurnal  period  (the 
minimums  at  t  ~  6.6,  11.8,  and  17.3  hrs)  induced  by  irregular  semidiurnal  (12.42  hr)  barotropic  tide. 
The  current  direction  tfft)  exhibits  complete  360°  rotation  during  approximately  12  hours  ( (f)  =  123°,  for 
example,  at  t  =  0.4  and  12.8  hr),  which  was  also  affected  by  higher  harmonics  (for  example,  <f)  =  270° 
at  t  =  6.8  and  17.3  hr  or  (j)=  180°  at  t  =  1.4  and  16.1  hr). 

The  MSS  profiles  £mc(C)  obtained  in  this  rotating  tidal  flow  were  compared  with  the 
dissipation  rate  profiles  ewl  (<f)  calculated  using  (i)  the  law  of  the  wall  and  (ii)  the  estimates  of  friction 
velocity  u,  deduced  from  logarithmic  approximation  of  the  ADCP  velocity  profiles  very  near  to  the 
seafloor  (£"  =  0.04  -  0.4  mab)  while  ensuring  the  skin-layer  accuracy  for  it, .  The  sensitivity  of  u,  to 
various  time  averaging  periods  was  also  tested;  30  min  averaging  provided  the  highest  coefficient  of 
determination  for  the  log-layer  fits,  while  ensuring  reasonable  stationarity  of  tidal  flow  over  the 
averaging  period.  During  the  tidal  cycle,  ut  varied  in  the  range  (l  -  7)x  10  m/s.  In  general,  law  of 

the  wall  dissipation  profiles  befittingly  intersected  with  the  MSS  profiles  at  2  -  5  mab,  suggesting 
the  applicability  of  the  law  of  the  wall  in  shallow  tidal-affected  waters  at  the  time  scales  of  about  30 
min  (see  examples  in  Fig.  4). 


Fig.  4.  The  MMS  dissipation  profiles  smc  (<f )  and  the  law-of-the-wall  profiles  sb  (tf )  for  the  across-slope 
eastern  (a)  and  along-slope  southern  (b)  flow  (rav,  =  30  and  10  min). 

When  turbulence  in  the  bottom  boundary  layer  is  influenced,  in  addition  to  local  friction,  by 
other  external  sources  (lateral  advection,  convection,  internal  wave  generation),  the  level  of  the 
dissipation  can  be  substantially  higher  compared  to  swl  (<f ) .  The  local  bathymetry  also  appears  to  play 

a  role  in  the  enhancement  of  smc  (tf )  compared  to  swl  {tf )  in  the  lower  7  -  10  m  of  the  water  column 
during  the  upslope  phase  of  the  westerly  directed  tidal  current  in  the  CDW  region  (Fig.  5). 
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Fig.  5.  The  MMS  dissipation  profiles  smc  (f )  and  the  law-of-the-wall  profiles  £b  (tf )  for  the  across-slope 
western  flow  influenced  by  a  turbulent  wake  behind  local  seamount.  (Tavr  =  30  and  10  min). 

It  is  hypothesized  that  the  turbulence  generated  at  the  summit  and  flanks  of  a  small  seamount  to 
the  east  of  the  observational  site  was  advected  onshore,  providing  an  increased  level  of  the  dissipation 
rate  in  the  BBL  compared  to  regular  background  tidal-induced  wall  turbulence. 

IMPACT/APPLICATION 

Our  research  program  has  been  enormously  strengthened  by  the  international  collaboration 
provided  by  the  grant.  A  joint  paper  with  the  Chinese  colleagues  has  been  published  in  2009  in  the 
Journal  of  Marine  Systems  and  another  paper  on  turbulence  intermittency  is  under  revision  in  the 
Journal  of  Geophysical  Research.  A  paper  on  the  near-bottom  turbulence  in  the  CDW  region  of  ECS 
is  ready  for  submission  to  the  journal  of  Continental  Shelf  Research.  A  complimentary  research  effort 
on  mixing  efficiency  in  stratified  flows  generated  a  paper  that  is  being  submitted  to  the  Journal  of 
Physical  Oceanography.  The  PI  Lozovatsky  visited  KORDI  in  August,  and  had  productive  discussions 
on  the  analysis  of  data  obtained  recently  by  the  Korean  collaborators.  He  will  be  visiting  China 
(Xiamen  University)  in  October  to  discuss  future  research  in  the  East  and  South  China  Seas  in 
collaboration  with  Chinese  oceanographers. 

TRANSITIONS 

None 

RELATED  PROJECTS 

The  Co-P.I.  Fernando  is  involved  in  another  ONR  funded  project  dealing  with  laboratory 
investigations  of  submarine  wakes  in  stratified  fluid  funded  by  the  ONR  turbulence  program. 
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